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Detecting Lagrangian fronts with favourable fishery conditions 
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Lagrangian fronts in the ocean delineate boundaries be- 
tween surface waters with different Lagrangian properties. 
They can be accurately detected in a given velocity field by 
computing synoptic maps of the drift of synthetic tracers, 
their Lyapunov exponents, and other Lagrangian indicators. 
Using Russian ship's catch and location data for a number 
of commercial fishing seasons in the region of the northwest 
Pacific with one of the richest fishery in the world, it is 
shown that the saury fishing grounds with maximal catches 
are located mainly along those Lagrangian fronts where pro- 
ductive cold waters of the Oyashio Current, warmer waters 
of the southern branch of the Soya Current, and waters of 
warm-core Kuroshio rings converge. Computation of those 
fronts with the altimetric geostrophic velocity fields both 
in the years with the First and Second Oyashio Intrusions 
shows that in spite of different oceanographic conditions in 
both the cases the front locations may serve good indicators 
of potential fishing grounds. 



1. Introduction 

In the northwest Pacific, the cold Oyashio Current flows 
out of the Arctic along the Kamchatka Peninsula and the 
Kuril Islands and converges with the warmer Kuroshio Cur- 
rent off the eastern shore of Japan. This frontal zone is 
known to be one of the richest fishery in the world due to 
the large nutrient content in the Oyashio waters and high 
tides in some areas therein. Some people attempted to iden- 
tify the favourable oceanographic conditions for catching 
pelagic fishes. By the common opinion, good fishing areas 
are often found at the boundaries of warm and cold cur- 
rents and around warm-core rings where the energy of the 
physical system is transfered in some way to biological pro- 
cesses. This strong physical-biological interaction provides 
favourable conditions for marine organisms. Surface conver- 
gent fronts of considerable physical and biological activity 
occur in zones where different water types impinge. Their 
admixture locally simulates phytoplankton photosynthesis 
and thus sustains higher phyto- and zooplankton concentra- 
tions, the main food of pelagic fishes. SST gradients have 
long been the main indicators used to find places in the 
ocean with rich marine resources. 

The Kuroshio- Oyashio fronts have long been recognised 
by Japanese fishermen to attract squids, fishes and mammals 
(for the earlier studies see Ref. [Uda, 1938]). It was found 
that fishing grounds depend on the location of the Oyashio 
fronts and vary from year to year [Fukushima, 1979; Saitoh 
et al., 1986; Sugimoto and Tameishi, 1992; Yasuda and Kita- 
gawa, 1996]. The fishing grounds are formed nearshore if 
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there exists the First Oyashio Intrusion along the eastern 
coast of Hokkaido, Japan. In the other years, the fishing 
grounds are located offshore where the Second Oyashio In- 
trusion is formed due to the presence there a large warm- 
core Kuroshio ring. It has been shown that locations of the 
fishing grounds depend not only on instantaneous and local 
oceanographic conditions nearby the fishing grounds but on 
the conditions over the whole region. 

We study here the connection of Lagrangian fronts (LFs), 
delineating boundaries between surface waters with different 
Lagrangian properties, with fishing grounds with maximal 
catches. To be concrete, we focus on fishing grounds of 
Pacific saury (Cololabis saira), one of the most commercial 
pelagic fishes in the region. Pacific saury, in general, migrate 
seasonally from the south to the north. In winter and spring, 
spawning grounds are formed in the south, off the eastern 
coast of Honshu, Japan. In spring and summer, juvenile and 
young saury migrate northward to the Oyashio area. After 
feeding in those productive waters, adult saury migrate to 
the south in the late summer. Commercial fishing begins in 
August and ends in December. 

Based on the AVISO altimetric geostrophic velocity fields, 
we compute synoptic maps of zonal, meridional, and abso- 
lute drift of synthetic tracers, their Lyapunov exponents and 
zonal and meridional entrance maps. Those maps with the 
catch data from the Russian saury fishery overlaid allow to 
identify the LFs in the region with favourable fishing condi- 
tions in the years both with the First and Second Oyashio 
Intrusions. It is shown that the LFs may serve a new indi- 
cator for potential fishing grounds. 

2. Data and Methods 

The method we used is based on the Lagrangian ap- 
proach to study mixing and transport at the sea surface 
[Haller and Yuan, 2000; Boffetta et al., 2001; Mancho et al, 
2004; d'Ovidio et al, 2004; Shadden et al, 2005; Kirwan Jr., 
2006; Lehahn et al, 2007; Beron-Vera et al, 2008; Prants 
et al., 2011a, b] when one follows fluid particle trajectories 
in a velocity field calculated from altimetric measurements 
and/or obtained as an output of one of the ocean circu- 
lation model. The important notion in that approach is 
so-called Lagrangian coherent structures (LCS) [Haller and 
Yuan, 2000]. The well-known LCSs in the ocean are eddies, 
jet currents, and fronts that can be visible in Eulerian veloc- 
ity fileds and at satellite images of the SST and/or chrophill 
concentration. There are another, less known LCSs, that 
are not visible at snapshots but can be computed with a 
given velocity field by special methods. These LCSs, for 
example, the stable and unstable invariant manifolds of hy- 
perbolic trajectories in the ocean can be identified with the 
help of some Lagrangian indicators (for reviews see [Wig- 
gins, 2005; KosheV and Prants, 2006]). The finite-size and 
finite-time Lyapunov exponents (FTLE) are the most used 
of them. Remarkably, the Lagrangian methods allow to 
identify some small-scale structures that are below the alti- 
metric resolution. These structures, typical in the processes 
of chaotic advection [KosheV and Prants, 2006; Budyansky 
et al., 2004], have been captured by computing the FTLE, 
absolute, zonal, and meridional drifts and other Lagrangian 
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quantities. The impact of LCSs on biological organisms has 
been recently studied in Ref . [ Tew Kai et al. , 2009] . By com- 
paring the seabird satellite positions with computed LCSs 
locations, it was found that a top marine predator, the Great 
Frigatebird, was able to track the LCSs in the Mozambique 
Channel identified with the help of the finite-size Lyapunov 
exponent. 

Geostrophic velocities were obtained from the AVISO 
database (http://www.aviso.oceanobs.com). The data is 
grided on a 1/3° x 1/3° Mercator grid. Bicubical spatial in- 
terpolation and third order Lagrangian polinomials in time 
have been used to provide accurate numerical results. La- 
grangian trajectories have been computed by integrating 
the advection equations with a fourth-order Runge-Kutta 
scheme with a fixed time step of 0.001th part of a day. 

The LCSs are determined mainly by the large-scale ad- 
vection field, which is appropriately captured by altimetry. 
It is expected that computation of particle's trajectories sta- 
tistically is not especially sensitive to imperfections of the 
velocity field caused by the interpolation and measurement 
imperfections. The comparison of the LCSs computed with 
altimetric velocity fields in numerous paper (see, for ex- 
ample, [Abraham and Bowen, 2002; Olascoaga et al., 2006; 
Beron-Vera et al, 2008; d'Ovidio et al, 2009; Huhn et al, 
2012]) with independent satellite and in situ measurements 
of thermal and other fronts in different seas and oceans has 
been shown a good correspondence. 

The FTLE is useful in detecting stable and unstable man- 
ifolds in irregular velocity fields. It is computed here by the 
method of the singular- value decomposition of an evolution 
matrix for the linearized advection equations introduced in 
[Prants et al, 2011a]. Of particular importance in detect- 
ing LFs is another Lagrangian indicator, the drift, D, that is 
simply the distance between the final and initial positions of 
advected particles [Prants et al, 2011b]. This quantity along 
with the zonal and meridional drifts has been shown to be 
useful in quantifying transport of radionuclides in the North- 
ern Pacific after the accident at the Fukushima atomic plant 
station [Prants et al, 2011b] and transport of the Madagas- 
car plankton bloom [Huhn et al., 2012]. 

The SST data (http://oceancolor.gsfc.nasa.gov) were 
used to illustrate oceanographic conditions in the cases of 
the First and Second Oyashio Intrusions. The data on fish- 
ing positions in latitude and longitude and daily catches 
were obtained from the database of the Federal Agency for 
Fishery of the Russian Federation. 

3. Results and Discussion 

3.1. Identification of the Lagrangian fronts favourable 
for saury fishing in the season with the First Oyashio 
Intrusion 

We restrict our analysis in this paper by the region to the 
east off Hokkaido and southern Kuril Islands coasts where 
the daily saury catch data from the Russian fishery were 
collected for 2001-2011. With the aim to detect the LFs 
separating waters of different origin and histories, we dis- 
tribute a large number of syntethic particles over that region 
and integrate advection equations backward in time for two 
weeks. We have computed different Lagrangian indicators, 
but only three of them will be used in the paper to visualize 
the LFs. 

In some years there exists the First Oyashio Intrusion 
along the eastern coast of Hokkaido island which is visible 
on the SST image in Fig. la averaged for 23-25 September 
2002. It is a combined map with the velocity field, con- 
tour lines of the FTLE, A, and locations of "instantaneous" 



hyperbolic (crosses) and elliptic (triangles) fixed points im- 
posed. Blue and red triangles mean cyclonic and anticy- 
clonic rotations, respectively. Down, \j , and up, A, tri- 
angles mean that the area of a tracer patch in that place 
decreases or increases, respectively. The data of saury catch 
for 23-25 September 2002 overlaid allow to conclude that 
the fishery grounds have been concentrated mainly in the 
waters of the First Oyashio Intrusion (the radius of the cir- 
cles in the figure is proportional to the catch in tons per a 
given ship). 

The FTLE is the finite-time average of the maximal sep- 
aration rate for a pair of neighbouring advected particles 
that can be computed by different methods. We used the 
method proposed in Ref. [Prants et al., 2011a] to compute 

A(r(i)) = ±W(G(t)), (1) 

where r is an integration time, a(G(t)) the largest singular 
value of the evolution matrix for linearized advection equa- 
tions. Scalar field of A is Eulerian but the very quantity is a 
Lagrangian one that measures an integrated separation be- 
tween trajectories. The ridges (curves of the local maxima) 
of that field delineate stable manifolds of the most influen- 
tial hyperbolic trajectories in a region when integrating ad- 
vection equations forward in time and unstable ones when 
integrating them backward in time. By definition stable 
manifolds specify two direction along which fluid particles 
approach a hyperbolic point whereas the corresponding un- 
stable manifolds specify two other directions along which the 
particles move away from that point [Wiggins, 2005; KosheV 
and Prants, 2006]. 

The motion around elliptic points is stable whereas it is 
unstable nearby hyperbolic points. Each hyperbolic trajec- 
tory 7(t) posesses stable W s (7(£)) and unstable W u {^f{t)) 
manifolds which are material lines consisting of a set of 
points through which at time moment t pass trajectories 
asymptotical to j(t) at t — ► oo (W s ) and t — ► — oo (W u ) (see 
[Wiggins, 2005; KosheV and Prants, 2006]). Those mani- 
folds are organizing structures in the flow, because they act 
as boundaries to fluid transport, attract and repel fluid par- 
ticles not belonging to them, and partition the flow into 
regions with different types of motion. The stable and un- 
stable manifolds coinside, as a rule, with LFs when inte- 
grating advection equations forward and backward in time, 
respectively. We would like to stress the important role of 
LFs because, in difference from rather abstract geometric 
objects like W s and W u , the LFs are fronts of real physical 
quantities that can be measured directly. 

Computing backward in time the displacements on the 
Earth sphere with the radius R 

D = R arccos [sin yo sin?// + cosi/o cos?// cos(x/ — xo)], 

for a large number of particles from their initial to final po- 
sitions, we get the absolute, D, drift map that enables to 
visualize clearly the fronts with convergent water masses of 
different origin and histories. In Fig. lb we show the FTLE 
field on 24 September 2002 with the contours of the particle's 
drift D. The black ridges on the map delineate the unstable 
manifolds which are a special kind of the LFs. We see that 
the saury fishing grounds with maximal catches are located 
mainly nearby the ridge with the maximal drift gradient. To 
understand why the fishing grounds correlate namely with 
this LF we compute the zonal, D x , and meridional, D y , 
drift maps on the same day. The color on the zonal drift 
map with the contours of A imposed (Fig. lc) distinguishes 
the particles entering to the region shown through its west- 
ern (red) and eastern (green) boundaries whereas red and 
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green colors on the meridional drift map (Fig. Id) mean 
that particles enter to the region through its southern and 
northern boundaries, respectively. It is seen on both the 
drift maps that productive cold waters of the Oyashio Cur- 
rent and warmer waters of the southern branch of the Soya 



Current flowing through the straights between the southern 
Kuril Islands converge at the LFs with the maximal catches. 
Animation of the daily Lagrangian maps for August - De- 
cember 2002 with the fishery grounds overlaid is available 
at http://dynalab.poi.dvo.ru/data/GRL12/2002 
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Figure 1. The season with the First Oyashio Intrusion on 24 September 2002. (a) SST image, (b) FTLE map, (c) zonal 
and (d) meridional drift maps with locations of saury catches imposed. FTLE is in units of [days] -1 and D x , y is in degrees. 



3.2. Identification of the Lagrangian fronts favourable 
for saury fishing in the seasons with the Second 
Oyashio Intrusion 

The oceanographic situation in the region cardinally 
differs in the years with the Second Oyashio Intrusion 



when a large anticyclonic eddy, formed as a warm- 
core Kuroshio ring, approaches to the Hokkaido eastern 
coast and forces the Oyashio Current to shift to the east 
rounding the eddy. 
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Figure 2. The season with the Second Oyashio Intrusion on 17 October 2004. (a) SST image, (b) FTLE map, (c) zonal 
and (d) meridional drift maps with locations of maximal saury catches imposed. 



The SST image in Fig. 2a averaged for 16-18 Oc- 
tober 2004 demonstrates that the fishery grounds have 
been concentrated mainly at the convergence front be- 
tween the Oyashio and Soya waters and the periph- 
ery of the anticyclonic Kuroshio ring with the center 
approximately at x = 147°30 / E, y = 42°30 / N. The 
FTLE map in Fig. 2b on 17 October 2004 demon- 
strates that the fishery grounds are located along the 
three main unstable manifolds in the region which 
are the LFs. The zonal and meridional drift maps 
on the same date in Figs. 2c and d show clearly 
that the fishery grounds with maximal catches are 



located at the LFs where productive cold waters of 
the Oyashio Current, warmer waters of the southern 
branch of the Soya Current and warm and salt wa- 
ters of the Kuroshio ring converge. Animation of 
the daily Lagrangian maps for August - December 
2004 with the fishery grounds overlaid is available at 
http://dynalab.poi.dvo.ru/data/GRL12/2004 
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In the autumn of the next year the oceanographic 
situation in the region has been characterized again by 
the Second Oyashio Intrusion but with the current sys- 
tem different from that in 2004 (see the STT image in 
Fig. 3a). In the FTLE map (Fig. 3b) on 17 Octo- 
ber 2005 we see the anticyclonic Kuroshio ring, similar 
to the ring in 2004, with the center approximately at 
xq = 148°E, 2/0 = 42°30 / N forming now a dipole pair 
with another anticyclonic warm-core eddy with the cen- 
ter approximately at x = 145°30'E, y = 40°30'N. 
Moreover, there is another dipole with two cold-core 
cyclonic eddies to the east of the first dipole. That 
four-eddy current system is clearly delineated by the 



black ridges of unstable manifolds in Fig. 3b. The 
zonal and meridional drift maps in Figs. 3c and d 
visualize the LFs of the current system and the LF 
separating waters of the Oyashio and Soya Currents. 
Figure 3 shows that the fishery grounds with maxi- 
mal catches are located now not only at the LF be- 
tween the Oyashio Current and the periphery of the 
Kuroshio ring, as in Fig. 2, but at the LFs around 
the four-eddy current system as well. Animation of 
the daily Lagrangian maps for August - December 
2005 with the fishery grounds overlaid is available at 
http://dynalab.poi.dvo.ru/data/GRL12/2005 
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Figure 3. The same as in Fig. 2 but on 17 October 2005. 



4. Conclusions 

We introduce the notion of a Lagrangian front (LF), 
the zone where surface waters of different origin and 
history converge. Based on altimetric velocities we in- 
tegrate advection equations for a large number of syn- 
thetic tracers backward in time and compute the FTLE 
maps and zonal and meridional drift maps in the re- 
gion to the east off the Hokkaido and southern Kuril 
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Islands coasts, one of the richest fishery place in the 
world. The data on fishing positions and daily catches 
of the Russian ships are imposed on the SST, FTLE 
and drift daily maps. It is shown that the saury fish- 
ing grounds with maximal catches are located mainly 
along those LFs where productive cold waters of the 
Oyashio Current, warmer waters of the southern branch 
of the Soya Current, and waters of warm-core Kuroshio 
rings converge. We have found some correlations be- 
tween locations of the thermal fronts with saury fishing 
grounds in the region but have not found any correla- 
tions with chlorophyll fronts. Thus, it is shown that LF 
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locations may serve good indicators of potential fishing 
grounds in rather different oceanographic conditions. 
The method proposed seems to be quite general and 
may be applied to forecast potential fishing grounds for 
the other pelagic fishes in different regions of the World 
Ocean. 
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